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Abstract 

The quantum dynamics of an electron in a uniform magnetic field is stud- 
ied for geometries corresponding to integrable cases. We obtain the uniform 
asymptotic approximation of the WKB energies and wavefunctions for the 
semi-infinite plane and the disc. These analytical solutions are shown to be in 
excellent agreement with the numerical results obtained from the Schrodinger 
equations even for the lowest energy states. The classically exact notions 
of bulk and edge states are followed to their semi-classical limit, when the 
uniform approximation provides the connection between bulk and edge. 
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I. INTRODUCTION. 



The aim of this work is to present some analytical methods to obtain the energy spec- 
trum and the eigenfunctions of non-interacting electrons constrained to a finite domain with 
boundaries and submitted to a uniform magnetic field. 

This problem is relevant to various situations in condensed matter physics. In the low 
magnetic field regime, defined by the condition $ -C $05 where $ is the magnetic flux 
through the system and $0 = ~ is the quantum flux, recent experiments performed on small 
metallic systems did show how important is the effect of the boundaries It determines 
the nature of the zero-field classical motion being either integrable or chaotic. The magnetic 
susceptibility has been shown using numerical and semiclassical methods to be reduced and 
to present large fluctuations in the chaotic case, whereas it is larger and varies slowly with 
the field in the integrable case 0,^. In the opposite limit $ ^ $0 of high fields, we are in the 
so-called integer quantum Hall effect regime (IQHE), where the edge states associated with 
the boundary play a prominent role In this work, we shall concentrate on the problem 
of non-interacting electrons in a semi-infinite plane and in a disc. 

The classical dynamics allows for a natural distinction between bulk and edge states. 
A first semiclassical method is based on the Einstein-Brillouin-Keller (EBK) quantization 
rules and preserves this bulk and edge states splitting, by giving different quantization 
rules for each of them. This approximation is further improved for the semi-infinite plane 
by constructing the asymptotically matched WKB function and then finding its zeros cor- 
responding to the energy levels. This procedure removes the gap between the bulk and edge 
energies, resulting in a very good approximation for the exact spectrum. The calculation is 
not new 0, however we present it in detail, as it serves a starting point for the WKB approx- 
imation for the disc. The spectrum of the disc is found using the more general comparison 
equation method (for a general reference, see 10), also called the Miller-Good method 
to map the problem onto the semi-infinite plane's one. The obtained semiclassical formulas 
are valid for any strength of the magnetic field, thus providing us a way to study more 
particularly both the IQHE and the low field regimes. 

II. THE WKB SPECTRUM OF THE SEMI-INFINITE PLANE 

In this section the spectrum of an electron in the semi-infinite plane in a magnetic field is 
approximated first using the EBK quantization rules and then building the matched WKB 
wavefunction. To introduce semiclassical language, we begin by considering the classical 
dynamics. 

A. The classical dynamics. 

We consider a spinless particle of charge — e (e > 0) and mass m constrained to move in 
the semi-infinite plane. A uniform magnetic field B is applied perpendicular to the plane. 
Cartesian coordinates are defined such that the x axis is perpendicular to the boundary and 
the motion is confined to positive values of x. It is convenient to consider the boundary 
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having a finite length L and therefore we impose periodic boundary conditions in the y 
direction so that the particle moves on a semi-infinite cylinder (Fig. Q). 
In the Landau gauge A = (0 , Bx), the Hamiltonian of the particle is: 

H=±^{p,'HP,^'-B.r), (1) 

and the momentum is p = {mi ,my — -Bx). The total energy E and the y-component 
Py of the momentum are constants of motion, therefore the problem is integrable. In the 
four-dimensional phase space of the Cartesian coordinates and the corresponding momenta, 
each family of classical trajectories are winding on an invariant torus defined by the two 
constants of motion. 

The ensemble of trajectories splits naturally into two families: those that do not touch 
the boundary (bulk trajectories), and others (edge trajectories). Bulk particles of energy E 
go counterclockwise in circles of radius = \J2E /muj'^ (where u = eB/mc is the cyclotron 

frequency) with their center farther than from the edge and momenta Py < —\j2raE. Edge 
particles have momenta Py > —\/2mE and undergo specular reflections on the boundary 
before closing a circle so that the cyclotron orbit center begins drifting along the edge 
(Fig. |l|). Increasing py at fixed energy, the bulk tori in phase space are transformed at 
Py = —V2mE in a discontinuous way into edge tori. 



B. The EBK quantization. 

We consider now the quantum-mechanical version of the same problem. Dirichlet bound- 
ary conditions are imposed on the wavefunction \I'(a;,?/): 

^(0,y) = 0. (2) 

The application of the EBK quantization rules for integrable systems leads to the quan- 
tization of the two actions ly and along the y and x axis. The action ly corresponding to 
the motion along the y axis (parallel to the boundary) is 
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]-jpydy = ^ = Hyh, (3) 



with Uy G TZi. The motion along the x axis (perpendicular to the boundary) is different 
for the bulk and the edge states, therefore their quantization differs too. In particular the 
energy of the bulk states is found from 

1 /" 1 /" I 6 E 1 

Ix = — f Px dx = - / J2mE - {py + -Bxy dx = — = h{nx + -), (4) 
271 J vrJV c uj 2 



where Ux € IN. For the edge states the EBK condition is 

Ix = - J J '2mE - {py + -Bx)^ dx = — Urccos - ^o\/l - ^oj 



= + h, (5) 
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where 



Py _ '^^riyL 



mure y/2mE 

For the bulk trajectories (^o < ^1); the Maslov index is | and we obtain degenerate Landau 
levels E = huj{nx + |) which correspond to states that do not feel the boundary. For the 
edge trajectories (|^o| < 1), the integration range is restricted because of the reflection on 
the boundary and the Maslov index is |, since there is one turning point (Maslov index 
|) and one reflection (Maslov index | associated to a change of sign of the wavefunction) . 
The energies are implicit solutions of (^, they are non degenerate and bounded below by 
huj{nx + |) for each n^. Note that E ^ oo when C,o ~^ 1- There is a singularity in the EBK 
spectrum, separating bulk and edge energies. 



C. Matching the WKB wavefunctions. 

Once the motion in the y-direction is integrated, the Schrodinger equation together 
with the boundary condition (0) reduces to a one-dimensional Sturm-Liouville problem. A 
systematic WKB analysis is well-developed for those kind of problems and improves the 
EBK quantization. 

We introduce the dimensionless variable x = \/2x/Ib, and set fo = V^Py/mculB and 



e = hu!/2E, where Ib = yhc/eB is the magnetic length. Using (|ID, the Schrodinger equation 
for the one- dimensional wavefunction ^(x) = e~*^^^/'* \E'(x, reads: 

% + ^ + ^o)' ) V^(5) = 0. (7) 



A subsequent change of variable £ = (x + Xq) — 1 = — + fo — 1 gives: 

eV"(0 - (e' + 20 /(O = 0, (8) 

where /(O = ^{x). Eq. ([7p is a Weber equation [§,|n|], and its solutions (vanishing at 
infinity) are: (p{x) = cDe__i{x + xq), where c is a constant and D^{u) a parabohc cylinder 
function. The Dirichlet boundary condition (|^) reads in the new variables: 

^(0) = /(^o - 1) = 0. (9) 

Since the solutions / of Eq. (H) depend only on e, without this condition the energies would 
be proportional to huj (Landau levels). However this condition makes the rescaled energies 
^ to depend on the (energy-dependent) parameter ^o- They are therefore given by an 
implicit equation of the form: 

where is given by (|^), and the label refers to different 'energy bands' (in the infinite L 
limit, hn^{C,o) gives the appropriate energy bands by continuously varying ,^o)- 
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Our purpose in this section is to use a semiclassical approximation in order to find explicit 
analytical expressions for the energies. As noted by Isihara and Ebina 0], for small e (i.e., 
large energies), (|^) is a standard example of equation where the WKB method is applicable. 
It has two turning points at = and ^2 = ~2. Sufficiently far from them, the WKB 
function is given by the following asymptotic expression 0: 

fwKB{0 = ie + 20~'^ (c,e-^l'''^' + c,e^l'''^'], (11) 



where Ci and C2 are arbitrary constants. In the vicinity of the turning points this approxi- 
mation breaks down. The potential (^^ + 2^) is then linearized in Eq. (^ and the resulting 
Airy equation can be solved exactly. One obtains the approximate solution fwKsiO 
matching the corresponding Airy solutions for each turning point with the WKB solutions 
(p!T|) valid to the right, between, and to the left of them. 

The function fwRsiO consists of five branches, related to the five overlapping intervals 
in which the WKB and Airy approximations are valid. We start from large ^'s. 



i. The first branch has a domain defined by ^ ea, where the expression ([Tl|) holds 
[^]. The vanishing of the wavefunction at infinity fixes C2 = 0, and we set arbitrarily 
Ci = 1. Thus, the first branch of the WKB function is: 

/aB(0 = (e^ + 20-^e-^/o^^*^. (12) 

It is always nonzero, so there are no solutions of Eq. (||) for ^0 — 1 . 

ii. The second branch approximates the exact solution near the first turning point, i.e. 
for 1^1 <^ 1, where + 2^ can be linearized and replaced by 2^ in equation (||). It 
consists therefore of a linear combination of two independent Airy functions Ai(t) and 
Bi(t) with t = ^57j^. When ^ is such that ei ^ ^ ^ et, t is large so the asymptotic 
approximations of Airy functions can be used, and we can approximate the exponential 
in (|12D by exp(— ^2). We match for such ^'s the solution of the linearized equation 

with fivKsiO ^^"^ determine the unknown constants of the linear combination. This 
gives: 

/^^L(0 = |^Ai(40- (13) 

(2e)6 63 

The energies are found by imposing the Dirichlet condition on this solution: 

^ ' """^ ' 0<l-eo«l, (14) 



hcu 2^0-2' 



where the a^^ are the n^th zeros of the Airy function Ai{t). These energies are 
associated with the 'whispering gallery' orbits (see Fig. |I]) which are concentrated 
in a very narrow region near the boundary. We notice that the energies (|1^) can be 



obtained from a 'generalized EBK rule' by letting be such that + 1 = ^(— a 
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in (g), instead of an integer. Note that there are no energies for ^0 > since the real 
zeros of Ai(t) are negative. 
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iii. Between the two turning points and ^2, the expression (pi]) holds again and to 

(2) 

determine the constants ci and C2 we do the matching with fwRsiO t^e interval 

2 2 2 2 

—65 <^ ^ <^ —63. This gives for <^ —63 and ,^ + 2^63 : 







fwKBiO = 2(-e' - 20"^ sin(^ J dt V-t2 - 2t + j). (15) 

The argument of the sine for ,^ = ^0 — 1 is vr ( ^ + ^ ), where is the classical action for 
the edge motion calculated in (^. Thus as expected we obtain the EBK quantization 
(^ for ,^0 between but not too close to 1 and —1. The comparison of the energies 
derived within the WKB approximation with the exact ones found numerically is 
shown in Fig. Ea. 



iv. This branch represents the function in the vicinity of the second turning point ^2 = 
—2. Repeating the scheme of matching with the previous third branch, we obtain for 
le + 2|«l: 

flv^KsiO = ^ f sin(f )Ai(-^(e + 2)) + cos(f )Bi(-^(e + 2)) | . (16) 

(26j6 V 

The energies correspond here either to edge trajectories which nearly complete full 
circles before being reflected or to bulk trajectories very close to the boundary (fig. 
|I|). In Fig. ^ we see that, for the lowest band {e ~ 1), the energies of the fourth 
branch, found by setting fwxBi^o ~ 1) = 0, become closer below ^0 — —0.5 to the 
exact energies than those of the third branch. 

V. The fifth branch fwRsiO is derived in a similar way. For — ^ — 2 ^ 6^: 

/&B(0 = (f + 20- (sin(^)e-r-v^ + 2cos(^)e^r-v^). (17) 

The equation to be solved for the energies is: 

1 71 



tan — = -e^i-^"'^osH~io)~(o y/^) ^ _ 1 > . (ig) 

When 1 is large, the right-hand side of this equation is a large number so that the 
argument of the tangent in the left-hand side is to a good approximation vr time an 
half-integer. Thus for |^o| sufficiently large the energies are very close to the Landau 
levels given by (^, i.e. 1/26 = n^^ + 1/2. The solutions of the energy equations for the 
fourth and the fifth branches and the exact energies are shown in Fig. EB. 



The graphs in Fig. 2a, and 2c show a very good agreement between the exact calcu- 
lation and the matched WKB approximation even for the lowest energies. The distinction 
between the bulk and the edge states disappears and the spectrum rises smoothly from ^ 
(from above) to infinitely large edge energies. 
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III. SEMICLASSICAL SPECTRUM OF THE DISC IN MAGNETIC FIELD 



We study here the similar problem as that of the previous section but for a disc of 
radius R and a magnetic field applied in the 2;-direction perpendicular to the disc. For this 
geometry, a convenient choice for A is given by the radial gauge A = ]^Br eg, where (r, 9) are 
polar coordinates and e^, the corresponding unit vectors. The two momenta canonically 
conjugated to r and 6 are pr = mr and pe = mr'^{6 — the 2;-component of the angular 
momentum. In the symmetric gauge the Hamiltonian of the particle reads: 

Since pg and the energy are conserved by reflexion on the boundary, we again deal with an 
integrable system. 

We can as before distinguish between edge and bulk trajectories: for —muj^^R'^ + Rrc) < 
Pe < —muj{^R^ — Rr c) the particle bounces and drifts along the boundary (edge trajectories); 
for —mujiJ^R^ — Rr^) < < it follows a closed-circle bulk trajectory; for pg > ^ or 
Po < —muj{\R^ + Rrc), the motion is forbidden or outside the disc. Note that positive values 
of angular momenta are associated with bulk trajectories whose cyclotron orbit centers are 
closer than 2rc to the origin. Let us emphasize that if R is of the order of or less than the 
magnetic length Is-, the distinction between edge and bulk motion becomes meaningless from 
the quantum-mechanical point of view, all cyclotron radii being greater than Ib-, and then 
also than R {huj/2 is a lower bound for the energy). For such low fields or small systems all 
classical trajectories of interest are of the edge type. The case of weak and zero fields has 
been extensively studied before PJT(]|. We will return to the weak field limit later to show 



the limitation of our methods in that regime. 

We now formulate the quantum-mechanical problem. Since the system is invariant by 
rotation along {Oz), we choose wavefunctions \l/(r, 6') = y9(r)e*'^. We introduce the dimen- 
sionless variables f = — , e = ^ and L = — ^ = le. The wavefunction (p(C) then satisfies 
the Schrodinger radial equation: 

We impose the boundary conditions lim^^o viO < </?(—) = 0. Note that 

equation ( pO]) with these boundary conditions has an exact solution given by ip{C,) = 

c^l'le~fc iFi (i±|tl!i _ 1 _|_ |/|- where c is a constant and iFi (a, c; m) the confluent 
hypergeometric function [11|. The semiclassical approximation of Eq. (pO|) becomes better 
as the magnetic field is decreased. Indeed, the energies E tend to their finite zero-field values 
when S — »• and therefore e = hu/E approaches zero. 

As a result of the singularity at the origin due to the centrifugal potential ^ in equation 
(pop, the WKB method fails when applied directly. To overcome this difficulty we follow 
Langer [|T2| and make the change of variable x = — In(^), where a = 2{l—lz). The resulting 



equation for the function f{x) = (p{^) is: 



;^)'0 = g(x)/(x), Qix)=e-'^-2e-^ + {'^r. (21) 
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This equation can be solved using a WKB analysis similar to that of section 2. The two 



turning points are x± = — ln(l =p c), where c = y 1 — (^)^- Note that in two dimensions 
(the case we consider here) there is no need for a Langer substitution as in the case of the 
radial equation for the hydrogen atom [|13|. For > |, there are no real turning points, 



thus the WKB function is like in (0) an exponentially increasing function when x — —oo; 
therefore there are no solutions for such /^'s, and we restrict our discussion to the classically 
allowed region h < 

However we shall prefer a uniform approximation approach to the WKB analysis, be- 
cause the two turning points x+ and x_ coalesce when — > | and Iz —oo (the value 
/z = I corresponds to a cyclotron trajectory centered at the origin; 1^ — oo is actually 
not relevant). Iz = is also a special case where one of the turning points (x-) goes to 
infinity. We make use here of the so-called comparison equation method, that gives correct 
approximations no matter how close are the turning points. This will enable us to map 
the problem of the disc to the previously solved semi-infinite plane case. In this method 
01, the approximate solution is written in terms of a function g{cr), solution of a simpler 
equation. We shall choose for this simpler 'comparison equation' the equation (^ we al- 
ready solved, in which we make for convenience the change of variable ^ = a/I a — 1 so that 
^ = (o"^ — j)g{(7), where e is a parameter to be determined later. The solution of (^I]) is 
obtained by multiplying g{a) by a slowly varying amplitude and allowing its argument a to 
depend weakly on x. In other words, we look for solutions of the form /(x) = A{x)g{a{x)). 
The equation obeyed by the mapping function a{x) is simplified by choosing A{x) = (^)~2 
and reads: 

/ '3' N 9 ^ / N / d'^ X 9 / 9 1 \ / dcr 1 (P da , i 

i^/Qi^) = {^f (a' - .) + (^)'^(^)-'. (22) 

If (^)~^ varies slowly enough or if the left-hand side is large enough, the second term in 
the right-hand side of equation (^) is negligible. This constitutes the approximation itself 
(for more details see 0). 

We define by analogy with (|^) the classical radial action variable: 



7r/r(xo) = <j> Pr dr = j dr = ^^^^'^ j y^— Q(a;) dx 



rayja(l~c) 



Xo 



— ^ ( - (e-^o - 1)2 + arccos( ) - Vl - arccos(^— e^^o - -) j (23) 



(if xq > x+, this action is purely imaginary and we use arccos(u) = —i arccosh(u), yl 
— z — 1 for all u > 1; if Xo < x_, we define Ir{xo) by integrating from xq to x_ so it 
is again purely imaginary and (^31) still holds with arccos(M) = iarccosh|u|, — u'^ = 
—i ^u^ — \ for u < — 1 ). 

For the above analysis to be meaningful, the mapping x ^ a has to be one-to-one, i.e. 
^ and ^ must be always non-zero. From (|2^ ) it follows that x± must map into This 
enables us to integrate our approximation of (|22|) , which takes the form: 

/r(xo) = 4( = o-(a;o) ) = -^(arccos^^o - ^o\/l - ^o^), (24) 

/vre 
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where Ix{Co) is the classical action (|^) for the semi-infinite plane. The energy ^ is obtained 
by letting y/^a{x^) = —1 in (^^. Taking 6{l) = for / < and 1 for / > 0, we obtain 

— = L-ie(l). (25) 

The mapping function takes the following forms for x —>■ ±oo and x —>■ x±: 

^/ia(x) ~ (e|/|x)2 if s ^ 1 and x ^ 



^ecrfa;) ~ — J — e~ 2 if a; ^ — 1 




v^a(x) T 1 ^ (y^(T^^) ' + 1 T c) if |e-^ - 1 ± c| « c. (26) 

The semiclassical wavefunction is then given (up to an arbitrary multiplicative constant) by: 



fwKB{x) 



''''^''^ ^ fwKBU=Vea{x)-l), (27) 



Qix) 

where fwKsiO is the WKB solution of equation (§). Moreover, since limj-^oo cr(a^) = +oo, 
it satisfies the same boundary conditions as in section 2, i.e. \im^^^ fwRsiO < °° 
fwKBi^o - 1) = 0, with: 

11- 

Xo = -ln( — ) = ln(— --— ); ^o = V^cr{xo), (28) 



c 



where = is the total magnetic fiux through the system in units of the fiux quantum. 

To find the set of discrete levels one has therefore to solve the implicit equation 

^^ = hnA^o) + l0{l), (29) 

where depends on e through (^8]). /in^(^o) is the band energy curve ([l^) for the semi- 
infinite plane, and also the band energy curve for the disc in the limit of infinitely large iV$ 
if Iz < 0. In fact, in this limit in dH) varies continuously when / runs over all negative 
integer values (a is continuous). We distinguish the following typical cases: 

i. A^<i> ^ 1, bulk states. 

For absolute values of the angular momentum |/| <^ and energies low enough for 



e = 1^ ^ 1 to hold, xq and are large negative numbers (see (^), (^Sf ) and (^ 



It follows from the analysis of the previous section that hn^{^o) is close to a Landau 
level fix + ^, i-e. levels with small angular momenta |/| <^ practically do not feel 
the boundary. For negative /, the level is close to the n^'th Landau level, for positive 
/ the level is close to the {fix + O'th Landau level, i.e. ^ ~ (n^ + |) + ^ , ""-x being a 
positive integer. One finds for 1^ the upper bound 1^ < 2n\\ (which gives back the 
classical condition for oo). The fact that there are no bulk levels with positive 

angular momenta I > Hx near the ra^^'th Landau level is seen clearly on the numerically 
calculated spectrum as a function of the magnetic field, by following the levels from 
-8 = (see Fig. |^ or P|). Each time the band index Ux is decreased by one, one level 
of positive / is removed. 
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ii. ^ 1, edge states. 

These are associated with negative values of /, |/| > A'"^. Using the results of section 
2, we see that there are no levels if > 1; i-e- if < 1 — c (since a is an increasing 
function). This corresponds to the lower classical bound for the angular momentum 
h > {lz)min, where {lz)min = -N<s.e - ^/2Nq>e. For ^ 1", the rescaled energy ^ = ^ 
becomes infinite. The high energy part of the edge spectrum in a given band is thus 
given by those Iz = le very close to (but higher than) {lz)min, i-e., by 

In other words, the energy E (for fixed R) varies like {\l\-N^f with the magnetic field 
for such |/| ^ (see Fig. H). Note that |/| may become infinite but this corresponds 
to Iz 0. Decreasing |/|, we go towards lower energies. In the limit ^ 1, A^^^e and 
1 = become large, and to calculate the lower energy part of the edge spectrum we 
approximate Xq ~ ln(]|^) and a ~ 2|/|e ~ ^. Finally, by decreasing further |/| below 
N<!,, the states transform gradually into Landau-like states located far apart from the 
boundary, with energy levels closer to the Landau levels. 

In the high field limit, so that edge states do correspond to large angular momenta 
with respect to h, the semiclassical analysis is expected to work well. This is expressed 
in equation ( pT]) by the fact that the small parameter is ^ ~ jfj- ~ For these states 
1/N<s> may be thought as an effective Planck constant. In Fig. ^ we show that the 
semiclassical bulk and edge energies are effectively very close to the exact energies for 
= 20. 

iii. Ncs, <^ 1: low field limit. 

When — s> 0, the exact energies of the states of angular momentum Ih tend to their 
zero-field values Eo = 2^^, where ji are the zeros of the Bessel function of order / 
(note that due to the time-reversal symmetry, states with opposite / are degenerate in 
energy at 5 = 0, i.e. ji = j-i). For small fields and fixed ^, e = |f = 0{uj) and Iz = le 
is small. By expanding the action Ir{xo) in (p3D up to second order in the field we find 

——Ir{xo) = 2^zq - 1 - 2arccos(2;o 2) + [bz — ) 

h\l\ ^Jzq - 1 / 

, ,x 3 / , 3 2 ..dz"^ N,^5z. 
+(.o-l)-((-l + --^)(--— ) 

1^2 19 8 \ 

+ ^(-^0 - 3 + - - 4) + 0{uj'6z^l (31) 

6L Zq Zq I 

where z = ^ = ^^^g^^ lim^_o z = zq, bz = z - zq and i + j > 3. Using (0), 



ID, (p5D and lim^^Q elr^xo) = 0, we obtain limt^^cCo = 1- Therefore we can use 



the 'generalized EBK rules' developed in ii. of section 2 to quantize the action /^(xo). 



Taking = in (pi]), we obtain that the square root energy ^/Pzq = J '^"^^l^° is equal 



to the first term of the uniform approximation of the zeros of the Bessel function of 



order /, namely it satisfies (see for ex. ||Tl[|): 
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= |/|v/i^- |/|arccos(zo"^) = ^("^nji (32) 

Our semiclassical analysis gives thus the exact zero-field energies up to higher terms 
in this expansion (i.e. with an error of the order of one percent for the lowest energies 
- including those with 1 = 0). This error decreases as a power law by increasing |/| or 



ji |jTT|. The EBK rules to quantized the classical action /^(xo) (substituting n{nx + |) 



in the right-hand side of Eq. (|32D) give already very good approximations to the exact 



energies [|T0| in the zero field case. We conclude that then the corrections due to proper 
treatment of the whispering gallery orbits are small. 

For small but non-zero u, the action Ir{xo) remains quantized through the same rules 
and is equal (in this approximation) to its zero field value (^2]). It immediately follows 
from (^) that Sz = + 0{uj'^). Replacing this value into the action we obtain the 
energies up to second order in uj: 



= jf + 21N, + M(i + ^) + 0(^3) (33) 



This result agrees with recent semiclassical perturbative calculations at low fields ||14 
However it does not coincide with the exact second-order perturbation theory result 
of Dingle |jl5|, where a different factor multiplying A''! in (|33| ) was obtained (the two 



expressions are equal in the limit of infinitely large energies) . This discrepancy should 
affect the magnetic susceptibility, which is a measure of the field dependence of the 



levels, for low Fermi energies. For sufficiently large energies (|33|) can be used to study 
the level crossings at low A^^'s. The values iV$ at which there are level crossings are 
given by the solutions of the second degree equation: 

^(1 - I) + 2N,{h - h) + jI - jI = 0, (34) 

where j/^ and ji^ are zeros of the Bessel functions of order respectively /i and I2. 
For arbitrary and I we can compute the WKB wavefunction using the results of the 



previous section and the formula (|2^ and (^). It is given by (rj is a small parameter): 
fwKBix) = Q{x)--* e^-^'-(^), if - e-^ + 1 - c » 

"""^ ^ ^ -c)j, if |e---l + c|«c 

1 + 77, -e""' + l + Or] 

^ / 

, sin( — ) Ai 




+ cos(^) Bi 



a\plc 




'2e' vMe(l + c) 

Q(a;)^3 |sin(^)e-^"''-(") + 2cos(^)e^"''(")|, if e"" - 1 - c> r/. (35) 
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Suppose that one introduce an Aharonov-Bohm flux line through the origin of the disc, 
and let us study the transport of charges induced by varying the flux $. The vector potential 
is changed by A — Aq + ^^e, where Aq is the vector potential created by the uniform 
magnetic field. The Schrodinger equation for this system is still of the form (^) but with / 
replaced by / + a, where a = $/$0; therefore all our results still apply modulo this substi- 
tution. The spectrum is identical for a's which differ by integer values (gauge invariance). 
Changing continuously the flux $ is equivalent to an electromotive force (e.m.f.) through 
the system. Consider a varying continuously from to 1. Then all eigenstates and energies 
of the system evolve adiabatically so that edge levels lower their energies and fall (at a = 1) 
into the level situated just below them. Bulk levels are nearly not affected except those with 
positive /, whose energies E = hu{nx + / + a + |) is increased by hu, the electron moving 
from one Landau level into another (see Fig. p. Introducing the Fermi energy Ep for an 
independent electron gas at zero temperature, we see that the effect of changing the flux 
by one unit is to remove electrons from the center of the disc and to transfer them to 
the boundary if Ep is between the n^j'th and the [rix + l)'th Landau level (as in Laughlin's 
argument |^). This conclusion is not specific to the disc geometry, since it is due to positive 
I states which are almost insensitive to the boundary when >> 1. This scenario is 
therefore generalizable to arbitrary billiards in the strong magnetic fields. 



IV. CONCLUSION 

We have presented a calculation of the energy spectrum of an electron in a perpendicular 
and uniform magnetic field for the half plane and the disc geometry. These situations are 
the simplest to describe features due to the presence of both perpendicular field and a 
boundary. In fact, explicit calculations can be made in both cases using the existence of a 
constant of motion resulting from the symmetry of the system (namely the momentum in 
the direction parallel to the boundary for the half plane, and the angular momentum for the 
disc). One can then reduce these two-dimensional problems into one-dimensional ones and 
use semiclassical methods to solve them. 

Because of the inherent splitting between bulk and edge classical trajectories, in the 
standard EBK approach one still have a clear separation between bulk and edge states. 
A more correct semiclassical approach removes this splitting by taking into account the 
coalescence of the caustics with the boundary which are due to edge orbits making nearly 
full circles before being reflected or to bulk orbits very close to the boundary. This removes 
the EBK gap between bulk and edge energies and results in a semiclassical spectrum which 
approximates surprisingly well the lowest levels, as far as we checked from our numerical 
calculations. Blaschke and Brack did independently reach similar conclusions using periodic 
orbits techniques in their recent work on the disc billiard [l^. In the quantum system 
under usual boundary conditions (like the Dirichlet case here considered), the notion of 
edge and bulk states is thus not anymore well-defined. To restore a precise meaning of such 
a dichotomy in the quantum case and describe for instance the transition from edge to bulk 
states in the context of the integer quantum Hall effect, one has to consider more general 
boundary conditions [jl^ pertaining to the family defined by Atiyah, Patodi and Singer []TB 



The extension of our results to general (non-integrable) billiards with smooth boundaries 
could be made or by using the more elaborate path-integration semiclassical methods, or. 
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for almost circular billiards, by expanding the wavefunction on the WKB disc wavefunction. 
Such billiards are interesting because of the existence in zero field of regimes of localization 
of the wavefunction in the angular momentum space |T^. Our approach of deriving the 



uniform approximation for the separable systems can be of relevance in connection with the 
recent work on uniform approximation of the diffraction contribution to the semiclassical 
spectral density for general magnetic field free billiards [pO[] . 
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FIG. 2b. The energy spectrum for = 0: circles are the zeros of the fourth branch, otherwise 
notations hke in Fig. 2a. 
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FIG. 2c. The energy spectrum for 
notations hke in Fig. 2b. 



0: triangles are the zeros of the fifth branch, otherwise 
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FIG. 4. The energy spectrum of the disc as function of the angular momentum (N^ = 20): 
circles are the zeros of the third branch, pluses - the zeros of the fourth branch and triangles are 
the zeros of the fifth branch. 
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